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Abstract 

We study the stability of the two-neutrino MSW solution of the solar 
neutrino problem, corresponding to solar transitions into sterile neu- 
trino, Ue —>■ Ug, with respect to changes of the total fluxes of and 
''Be neutrinos, <1>b and $bc- For any value of <I>Be from the interval 
0.7^^1 < $Be < l-3$Br (for $Be = 0.7$!^) the i^e Us MSW transi- 
tions provide at 95% C.L. a description of the existing solar neutrino data 
for 0.40 (0.39) ^ $b ^ 36 (40)$^^, and being the fluxes 
in the solar model of Bahcall-Pinsonneault from 1992. The corresponding 
allowed regions of values of the parameters Am^ and sin^ 26, characterizing 
the solar neutrino transitions, are derived. The physical implications of the 
found MSW Vg Vg solutions for the future solar neutrino experiments are 
considered as well. 
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1. INTRODUCTION 

Recent developments in the domain of the solar neutrino problem [1,2] (see also the 
review article [3]) brought two new important elements. First, the spectrum of solar model 
predictions [4-9] for the flux of *B neutrinos, ^b, now spans the entire range of estimated 
99.9% CL. uncertainties [2]. The central value of <I>b derived in the "low" flux model [9] is 
smaller approximately by the factors 2.0 and 2.3 than the central values in the "high" flux 
models of refs. [4] and [8], respectively []. Second, it was found that the upper limits on the 
^Be neutrino flux, $Be; which can be inferred from the existing solar neutrino data [1,12-14] 
(see also ref. [15]), are significantly lower [16] than the values predicted by the solar models. 
The predictions [4-9] for $bc vary only by ~25% (from $Be = 4.20 x 10^ z/g/cm^/sec in ref. 
[9] to $Bc = 5.18 X 10^ z/g/cm^/sec in ref. [8]). No plausible astrophysical and/or nuclear 
physics explanation of the indicated beryllium neutrino deficit has been proposed so far, 
which represents a major new aspect of the solar neutrino problem. 

Having in mind the spread in the current solar model predictions for $b and ^bo, "we 
investigate in the present article the stability of the MSW solution of the solar neutrino 
problem [17-19] with solar z/g transitions into sterile neutrino Ug (see, e.g., ref. [19]), with 
respect to variations of the values of total fluxes of ®B and ''Be neutrinos. Assuming that the 
^Be neutrino flux has a value in the interval O.T^b^ — "^Bc < 1-3$bc 5 where ^Be" is the flux 
predicted in the reference solar model of Bahcall - Pinsonneault [4], we determine the range 
of values of the ^B neutrino flux, for which the results of the solar neutrino experiments can 
be described in terms of two-neutrino MSW transitions of the solar neutrinos into a sterile 
neutrino, z/g —>■ Vg- Similar analyses for the MSW solution with solar z/g transitions into an 
active neutrino [18,19], Ve '^tJ-ir), and for the different possible vacuum oscillation solutions 

^Extensive discussions of the possible sources of uncertainties in the solar model predictions for 
<^B are given in refs. [2-5,8-11]. 
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[20,21,22] 0, were performed respectively in refs. [23,24] and [22] (see also [11,25]). We find 
that the MSW z/g Ug transitions provide at 95% CL. a solution of the solar neutrino 
problem for any value of from the interval O.T^b^ < $Be < l-3$if (for <l>Be = 0.7$b^) 
if $B lies within the remarkably wide interval 0.40 (0.39) ^ $b ^ 36 (40)^ f , 
being the neutrino flux in the reference model [4]. The corresponding allowed regions of 
values of the parameters Am^ and sin^ 29 , characterizing the transitions, are derived. The 
physical implications of the found "low" and "high" $b MSW Vs transition solutions 

for the future solar neutrino experiments are considered as well. 

We shall use the latest published data from the solar neutrino experiments [1,12-14] in 
the present analysis: 

R(Ar) = (2.55 ± 0.25) SNU, (1) 

= (2.89 ± 0.42) X lO'' cm-2sec-\ (2) 

RGALLEx(Ge) = (77.1 t?oA) SNU, (3) 

RsAGE(Ge) = (69 ± 13) SNU, (4) 

where R(Ar), and RGALLEx(Ge) and RsAGE(Ge), are respectively the average rates of ^'^Ar 
and ^^Ge production by solar neutrinos observed in the experiments of Davis et al. [1], 
and GALLEX [13] and SAGE [14], and $3^^ is the flux of '^B neutrinos measured by the 
Kamiokande collaborations [12]. In eqs. (1) - (4) the quoted errors represent the added in 
quadratures statistical (1 s.d.) and systematic errors. 

2. MSW ^ Vs TRANSITION SOLUTIONS: GENERAL PROPERTIES 

It is convenient to introduce the parameters 

fB-^>0, fBe-||^>0, (5) 



^More complete lists of references on the vacuum oscillation solution of the solar neutrino problem 
can be found in refs. [21,22]. 
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in terms of which we shall describe the possible deviations of the and ^Be neutrino fluxes, 
$B and $Be; from the values they have in the reference model [4]. The fluxes $b and $Be 
predicted in the models [4,5,8,9] correspond, respectively, to fB = 1.0; 0.78; 1.14; 0.50, and 
fBe = 1.0; 0.89; 1.06; 0.86. 

The Kamiokande data imposes at 99.73% (95%) C.L. the following hmits on the values 
fB can possibly have [22]: 



The lower hmit (6) holds [22,11] in the case of solar two-neutrino transitions or oscillations 
into active (^'//(r) or ^'/^(r)) or sterile {ug) neutrino, as well as for transitions or oscillations of Ug 
involving more than two neutrinos (active and/or sterile). The limit (6) is universal: it does 
not depend on the type of possible transitions, and on the specific mechanism responsible 
for them. 

Contrary to the lower limit (6), the upper limit (7) is valid only for two-neutrino solar 
Vg transitions or oscillations into an active neutrino or [22,11]. It does not apply, 
in particular, to transitions of the solar i/^ into sterile neutrino (s). 

Owing to the specific dependence of the MSW two-neutrino transition probability on 
the neutrino energy E it is possible to derive, using the Kamiokande or the Cl-Ar data, a 
somewhat constraining upper bound on fB in the case of — > MSW transitions. For 
sin^ 2^ ^ 4 X 10"^ [26] the relevant (averaged over the corresponding region of Ue production 
in the Sun, etc.) MSW Ug survival probability, Pmsw(E), satisfies [17] min Pmsw(E) = sin^ 9 
for certain interval of values of E when the transitions are adiabatic. Utilizing this property 
of Pmsw(E) and taking into account that the contribution due to the ^B neutrinos in the 
argon production rate R(Ar) in the Cl-Ar experiment cannot exceed the measured value of 
R(Ar), eq. (1), one finds: 



fB 0.30 (0.37), 



(6) 



fB ^ 4.8 (4.3). 



(7) 



MSW, I/e ^ 



s 



fB < 67 (61). 



(8) 
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This bound corresponds to min Pmsw(E) = sin^^ = 8 x 10~^, and represents the 99.73% 
(95%) C.L. hmit. For values of sin^ 9 different from the indicated one the suppression due 
to Pmsw(E) of the integral entering into the expression for the contribution of neutrinos 
in R(Ar) (see, e.g., [11]), as can be shown, is weaker and the upper bound on fe one gets 
is not an absolute upper bound. Actually, 8 x 10~^ is the minimal value of sin^ 9 for which 
one has Pmsw(E) = sin^ 9 for all neutrinos having energy E ^ 3.0 MeV. In the same way 
one can use the Kamiokande data to obtain an upper limit on fe, which, however, turns 
out to be less constraining: for sin^^ = 8 x 10~^ one obtains fe ^ 100 (85). Obviously, the 
maximal value of fe, allowed by the Cl-Ar data, will be smaller than that in eq. (8) (and 
may occur for a slightly different value of sin^ 9) if one takes into account the contributions 
due to the ^Be and the CNO neutrinos in R(Ar), as will be done in our analysis. 

The remarkable difference between the upper bounds on fe implied by the Kamiokande 
data in the cases of v^. — > ^^{t) and of Ue MSW transitions has the following origin. 

The maximal allowed values of fe correspond to maximal possible suppression of R(K) due 
to the MSW transitions, being inversely proportional to the relevant (effective) suppression 
factor (see, e.g., refs. [22,11]). Since the active neutrinos contribute to the Kamiokande 
signal, while the sterile neutrinos Ug do not, the signal due to the ^B neutrinos having energy 
E is suppressed in the cases of Ve ^ij.{t) and z/g transitions respectively by the 

probability factors [Pmsw(E) + 0.15 (1 — Pmsw(E))] and Pmsw(E), where the term with the 
coefficient 0.15 is due to the contribution to R(K) [22]. As a consequence of the specific 
dependence of Pmsw(E) on E/Am^ and sin^ 2^ the strongest possible suppression of R(K) 
due to the transitions v^. — > Ug is by the factor ~ (6 — 8) x 10~^, and it is by a much bigger 
factor ~ 0.15 when the transitions are into active neutrino i^,j.{t). 

In the present study of the stability of the MSW u^. — > Vg transition solution with respect 
to $B and $Be variations, we use the approach adopted in the similar studies of the MSW 
— ^ti{T) transition [23] and of the vacuum oscillation [22] solutions. Namely, the fluxes of 
the pp, pep and the CNO neutrinos (see, e.g., refs. [2]) are kept fixed and their values were 
taken from ref. [4]. The fluxes of the ^B and ^Be neutrinos, and correspondingly, fe and 
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fee, are treated as fixed parameters, which, however, are allowed to take any values within 
certain intervals. In the case of $Be the interval chosen corresponds to 

0.7 < fee < 1.3. (9) 

It is somewhat wider than the interval formed by the current solar model predictions: 0.86 - 
1.06. For $B values in the interval determined by the inequalities (6) and (8) were considered. 

The above approach is motivated by the fact that the contributions of the CNO neutrinos 
to the signals in all three types of detectors [1,12-14] are predicted to be relatively small 
[2,4-9], and that (apart from the CNO neutrinos) the spreads in the predictions for the 
fluxes <I>B and $Be are the largest. A more detailed discussion and further justiflcation of 
this approach are given in refs. [23,22] (see also [11]). We would like only to emphasize here 
that the aim of our study (as like of the studies of the MSW Ue ~^ ^ii(r) transition and of 
the vacuum oscillation solutions performed in refs. [23,22]) was, in particular, to determine 
the ranges of values of $b and $Be for which the possibihty of MSW Vg transitions of 
solar neutrinos cannot be excluded by the existing solar neutrino data. Certainly, values of 
$B corresponding to, e.g., fs = 30 cannot be expected at present to appear in any realistic 
solar model. 

In the absence of "unconventional" behavior (MSW transitions, vacuum oscillations, etc.) 
of solar neutrinos, the signals in the Cl-Ar and Ga-Ge experiments can be written in the 
following form within the above approach: 

R(Ar) = (6.20fB + LlTfee + 0.40cno + 0.23pep) SNU, (10) 

R(Ge) = (70.8pp + 3.1pep + 35.8fBe + 13.8fB + 7.9cno) SNU, (11) 

where 6.20fB SNU is the contribution in R(Ar) due to the neutrinos, etc. 

We have applied the "method in the present analysis. In computing the for a 
given pair of values of the parameters Am^ and sin^ 29 we have neglected the estimated 
uncertainties in the reference model predictions [4] for the solar neutrino fluxes since the 
ranges within which $b and $Be have been varied exceed by far the uncertainties. We did, 
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however, take into account the uncertainties in the detection cross-sections for the detectors 
[1,12-14]. 

A recent study [19] of the solar Ve —>■ MSW transition hypothesis, performed within 
the reference model [4] (fe = fee = 1) including the estimated uncertainties in the predictions 
for the pp, pep, ^Be, and CNO neutrino fluxes, has shown that at 95% CL. there exists 
a small mixing angle nonadiabatic solution with 

3.0 X 10"^ eV^ ^ Am^ ^ 8.1 x 10"^ eV^ (12a) 

3.3 X 10"^ ^ sin2 2^ ^ 1.3 X 10"^ (126) 

which provided a good quality of the fit of the data available by April 1995 Q (min = 1.65 
for 2 d.f.). A large mixing angle (adiabatic) solution (present in the case of z/g —>■ i^^(t) MSW 
transitions [18,19]) was shown [19] to be excluded at 98% CL. by the data. 

Allowing fe and fee to take values in the intervals (6), (8) and (9), respectively, we have 
searched for Ue — > Ug MSW transition solution in the region 10~^ eV^ < Am? < 10"^ eV^ 
and 10"'^ < sin^ 29 < 1. It was found that the values of fe for which the hypothesis of solar 
transitions provides at 95% CL. a description of the solar neutrino data form a 
wide interval, whose width depends somewhat on fee- We give below these 95% CL. allowed 
intervals for fs in the cases of fee = 0.7; 1.0; 1.3: 

0.39 ^ fe ^ 40, (13a) 

0.39 ^ fe ^ 37, (136) 

0.40 ^ fB ^ 36. (13c) 

The maximal values of fe in eqs. (13a) - (13c) are determined by the Cl-Ar data (1) and 
therefore exhibit a certain, although weak, dependence on fee- The value of fe is constrained 



fBe = 0.7, 
fBe = 1.0, 

fee = 1-3, 



^In July 1995 the GALLEX collaboration updated their results adding new data from 9 runs of 
measurements to the previously existing data from 30 runs [13]. 
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from below basically by the Kamiokande result (2). The corresponding allowed regions of 
values of the parameters Am^ and sin^ 29 he all within the intervals: 

3.0 X 10-^ eV^ ^ Am^ ^ 1.1 x 10"^ eV^ (14a) 

8.0 X 10-^ ^ sin^ 29 ^ 0.45. (14&) 

Note that the Am^— interval (14a) of all the solutions with 0.39 ^ fe ^ 40 is only slightly 
wider than the interval (12a) of the values of Am^ of the fe = fee = 1 solution. Our results 
are illustrated in Figs, la - Ic where 95% C.L. allowed regions in the Am^ — sin^ 29 plane 
for a number of values of fe from the interval 0.40 < fe < 35 and for fee = 0.7; 1.0; 1.3 are 
shown. The values of Am^ allowed for given fe and fee as well as the solution intervals (12a) 
and (14a) are determined primarily by the Ga-Ge data, while the solution values of sin^2^ 
are constrained by the Kamiokande and the Cl-Ar results. 

For any given allowed by the data value of fe < 3 the MSW u^, — > Ug transition solution 
is of the nonadiabatic (NA) type, i.e., it corresponds to nonadiabatic transitions of at least 
a large fraction (the higher energy) ®B neutrinos, while the transitions of the pp, ^Be, pep 
and CNO neutrinos are adiabatic as long as sin^ 2^ ^ 4 x 10~^; for smaller values of sin^ 29 
some of the pp and CNO neutrinos and the ^Be and pep neutrinos can undergo nonadiabatic 
transitions as well. This solution is the u^. — > Ug transition analog of the MSW u^. ^h{t) 
transition NA solution at fe < 3 discussed in detail in ref. [23]. The existence of the NA 
solution of interest, the location of the allowed region corresponding to a given value of fe, 
the movement and the change of the size of this region when the value of fe is changed, etc. 
can be explained in the same way as this was done for the case of the Ue ^ij.{t) transition 
NA solution in ref. [23]. 

As fe increases beyond the value ~ 2 the "dynamics" of the u^. Vg transition solution is 
the following. At fe ^ 3 an adiabatic (AD) solution appears in addition to the nonadiabatic 
one for 3.0 x 10~^ ^ sin^ 29 ^ 0.46 and for values of Am^ which he within the interval (14a) 
(see Figs, la - Ic). For 3 ^ fe ^ 20 the regions of the NA and of the AD solutions in 
the Am^ — sin^ 29 plane are disconnected. With the increase of fe the two regions approach 
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each other. They become continuously connected at fe ~ 22. As fs increases further the 
subregion of the NA solution diminishes. 

In contrast to the respective NA solutions at fe < 3, the ^^i(t) and the Vf. — > Vs 

adiabatic solutions differ considerably: the first occurs for [23] (see also the first article 
quoted in [18]) 1.0 ^ f b ^ 3.4, 0.15 ^ sin^ Id ^ 0.90 and 6.2 x 10"^ eV^ ^ Am^ ^ 10"^ eV^ 
while the second holds for f b ^ 3 and typically smaller values of sin^ 2Q and Am^, as Figs, 
la - Ic indicate. None of the allowed regions in the Am^ — sin^ 2Q plane corresponding to 
the two adiabatic solutions overlap (compare our Fig. la (lb) with Fig. lb (la) in ref. [23]). 

For given f b < 3 and fse from (9) the region of the NA Vf. — > Vg solution is shifted 
to smaller values of Am^ (on average by a factor of 1.2) with respect to the region of 
the NA Vf. — > Vjiij) solution corresponding to the same fs and fse values (compare Figs. 
6a (6c) with Fig. 6b (6d) in ref. [19] as well as our Fig. la (lb) with Fig. lb (la) 
in ref. [23]). This difference is due to the fact that in matter v^{j^ scatters (coherently) 
on electrons, neutrons and protons while Vg does not. As a consequence the probabilities 
Pmsw(^) PmswI^) associated with the v^. Vg and v^, — >• transitions depend 

on different density ratios [27], i.e., on (N° - N°/2)/N'^^^ and N^/N"^^^ respectively, where 
Ng and N° are the electron and neutron number densities in the point of Vf. production 
in the Sun, and = Am^ cos2^/(2E\/2Gf) is the resonance density. One has in the 
region of u^, production in the Sun [4-9]: N°/N° = (0.45 - 0.20). The difference between 
Pmsw(E) and P|;^sw(E) due to the finite value of the ratio N°/N° is insignificant when 
(N° -N0/2)/N''*^'' » l,tan2 2^, or if N'^JW^'' « 1. However, for the transitions of the 
^Be, pp, pep and CNO neutrinos neither of the above conditions is fulfilled, the difference 
between the two averaged MSW probabihties is not neghgible and leads to the shift in Am^ 
between the Ug —>■ Vg and v^, — > z/^(^) allowed regions. 

As we have mentioned earlier, for fixed fB and fBe the solution values of sin^ 2Q are 
determined primarily by the Cl-Ar and the Kamiokande data. In the case of the adiabatic 

~^ solution they can be easily understood quahtatively due to the fact that the 
neutrino flux is suppressed approximately by the (energy-independent) factor sin^ 9. As 
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a consequence, the allowed values of sin^^ for a given is, corresponding to an adiabatic 
solution, are constrained by the Kamiokande and the Cl-Ar data to lie in the interval (95% 
C.L.): 

0.37 ^ Ib sin^ e < 0.49, (15) 

where the upper limit follows from (1). For fe = 3, for instance, eq. (15) implies 
0.43 ^ sin^2^ < 0.55. The exact calculations for fee = 0.7 (1.0) give in this case: 
0.41 ^ sin2 26' ^ 0.46 (0.45) (see Figs, la - Ic). The maximal allowed value of sin^ 2^ 
is smaller than the one implied by eq. (15) because, in particular, the upper limit in (15) 
was obtained without taking into account the contributions of ^Be and (pep + CNO) neu- 
trinos in R(Ar). The upper bound on sin^2^ following from (15) is closer to the maximal 
solution value of sin^ 29 for larger values of fe, for which the indicated contributions in R(Ar) 
are smaller. In the case of, e.g., fe = 10 (20) from (15) one finds sin^ 29 < 0.19 (0.096), while 
our results (Figs, la - Ic) show that sm^29 ^ 0.18 (0.094). 

If 0.40 ^ fe < 1 the hypothesis of MSW Ue transitions of solar neutrinos provides a 
good quality of the fit of the data (1) - (4): for fee = 0.7; 1.3 we have min — 1-22; 1.27 (for 
2 d.f.) reached at fg = 1.0; 0.99, Am^ ^ 4.3 x 10"^ eV^ and sin^ 29 ^ 7.6 x 10"^. A somewhat 
better description of the data is achieved for larger values of fe: if we fix fee = 1.0 (0.7) 
and consider the interval 0.40 < Ib < 3.5, the minimal value of the x^— function occurs 
at fe ~ 3.4, Am^ ^ 4.3 (4.6) x 10-^ eV^ and sm^29 ^ 1.9 (1.8) x 10-^, and one has 
min = 0.85 (0.82). For the adiabatic solution with fe < 5 and, e.g., fee = 1-0, one finds 
min = 4.3 at Ib = 5.0, Am^ ^ 6.3 x 10"^ eV^ and sin^ 29 ^ 0.29. 

3. PHYSICAL IMPLICATIONS OF THE MSW ^ i^s SOLUTIONS 

We shall discuss next the physical implications of the MSW Vg transition solutions 
for the future solar neutrino experiments SNO [28], Super-Kamiokande [29], BOREXINO 
[30] and HELLAZ [31]. Consider first the "low" ^B neutrino flux solution, 0.4 ^ f b ^ 0.6, 
which holds for small values of sin^2^: 8.0 x 10^"^ ^ sin^26' ^ 4.5 x 10^^. For this solution 
the ®B i/g spectrum and the corresponding spectrum of the recoil e" from the reaction 
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u + e~ — > u + e~ due to neutrinos to be measured respectively in SNO and Super- 
Kamiokande experiments are predicted to be relatively weakly deformed (see Figs. 2a and 
2b). The deformation of the '^B z/g spectrum caused for given Am^ and sin^ 2^^ by the 
MSW z/g — * Us transitions practically coincides with the deformation produced by the MSW 
i^e '^Kt) transitions for the same values of the two parameters. The deformations of the 
recoil e~ spectrum generated by the two types of MSW transitions of '^B neutrinos differ in 
general. This difference can be quite substantial for relatively large values of sin^ 29, but 
diminishes with the decreasing of sin^ 29 and becomes hardly observable for sin^ 29 ^ Ax 10^^, 
i.e., for the "low" $b solution under discussion (compare our Figs. 2a and 2b respectively 
with Figs. 4c, 4d and 6c, 6d in ref. [19] and with Figs. 2 and 3 in ref. [23]). The ratio 
Rcc/NC of the rates of events due to the charged current (CC) and neutral current (NC) 
reactions Ve + D — > + p + p and v + D — > z/ + p + n by which the ^B neutrinos (with energy 
E > 6.44 MeV and E > 2.2 MeV, respectively) will be detected in the SNO experiment, is 
not sensitive to the transitions (or oscillations) into sterile neutrino 

Possible signature of the "low" $b solution under discussion could be i) a sufficiently 
strong suppression of the flux of 0.862 MeV ''Be neutrinos by a factor not exceeding, e.g., 0.15, 
or ii) a specific deformation of the spectrum of pp neutrinos having energy E ^ 0.22 MeV. 
The 0.862 MeV ^Be neutrinos are predicted to produce approximately 90% of the signal in 
the BOREXINO detector [30], while the HELLAZ detector [31] is envisaged to measure the 
total flux and the spectrum of pp neutrinos with E > 0.218 MeV. In both detectors the 
V — e~ elastic scattering reaction will be utilized for detection of solar neutrinos. The pp Ve 
spectrum will be reconstructed from the data the HELLAZ detector is conceived to provide 
about the recoil-electron spectrum in the e~ kinetic energy region Te > 0.1 MeV, and the 



^For the "low" $b MSW i^^i^r) sohition on has [23]: Rcc/NC - 0.75 (0.85)RSM^nc 

fe = 0.5 (0.4), where Rcc/NC ratio of the CC and NC event rates predicted in the absence 

of "unconventional' behavior of the '^B neutrinos. 
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measurement of the recoil-electron momentum direction 0. 

Since the BOREXINO detector is based on the reaction u + ^ u -\-e^ [30], the active 
neutrinos i^^{t) are going to contribute to the BOREXINO signal, while the sterile neutrinos 
Us will not. As a consequence, the event rate generated by the 0.862 MeV ''Be neutrinos 
in BOREXINO can be suppressed at most by the factor 0.21 if they undergo — > z/^(t-) 
transitions, and it can be reduced by a noticeably smaller factor when the transitions are 
of the Ue Us type. In the latter case the suppression factor coincides with the averaged 
MSW probability Pl^f^ = Pmsw(E = 0.862 MeV). 

If 0.4 ^ fe ^ 0.6 and we consider only the corresponding values of Am^ and sin^ 29 allowed 
(at 95% C.L.) by the data (1) - (4), 

Pmsw constrained to lie in the interval ~ (0.06 — 0.40), 
the minimal and maximal values of Plifw taking place at Am^ = 3.1 x 10~^ eV^ and 
Am2 ^ 8.0 X IQ-^ eV^, but at the same sm^2e = 4.5 x lO'^. One has Fli^ < 0.15 for 
5.3 X 10"^ eV^ ^ Am^ ^ 8.0 x 10"^ eV^. The lowest (highest) value of this Am^-interval 
depends on the value of sin^ 29: the one given corresponds to sin^ 29 = 10~^ and it changes 
from 5.3 (8.0) x 10"^ eV^ to 3.1 (7.4) x 10"^ eV^ when sin^ 29 increases from 10'^ to 4.5 x lO^^. 
For values of Am^ in the solution interval 3.0 x 10"^ eV^ ^ Am^ ^ 4.5 x 10~^ eV^ the 
MSW Ue — ^ Us transitions lead to a relatively large distortion of the spectrum of pp (z/g) 
neutrinos (recoil electrons) having energy E ^ 0.30 MeV (Tc ^ 0.16 MeV), which might 
be observable with the HELLAZ detector; for 3.0 x 10'^ eV^ ^ Am^ ^3.7 (4.0) x 10'^ eV^ the 
pp z/g (recoil-electron) spectrum will be strongly deformed in the entire neutrino (electron) 



^Prom the kinematics of the reaction u + e u + e it follows that for the initial neutrino (recoil- 
e" kinetic) energy E < 0.42 MeV (Te > 0.10 MeV) one has Te < 0.26 MeV (E > 0.218 MeV). This 
determines the energy intervals in which the two spectra are expected to be measured with HEL- 
LAZ: 0.218 MeV ^ E ^ 0.42 MeV and 0.10 MeV ^ Tc ^ 0.26 MeV. 

^For the same values of A.m? and sin^ 29 the pp neutrino spectrum will be deformed, but much 
weaker, by MSW Ve — *■ ^/^(r) transitions. 
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energy interval E = (0.218 - 0.42) MeV (Te = (0.10 - 0.26) MeV) relevant to the HELLAZ 
detector Q This is illustrated in Fig. 3a and Figs. 4a - 4c, where the pp z/g spectrum 
and the spectrum of the recoil electrons from the reaction u + e~ —>■ u + e~ induced by 
the pp neutrinos, deformed by the MSW z/g Ug transitions under discussion, are shown 
for Am^ = 3.0; 3.5; 4.0; 4.5; 5.0 x 10"*^ eV^ and respectively for sin^26' = 5 x 10"^ and 
sin^26' = 10"^, 5 x 10"^, 10~^, 2 x 10"^. As our calculations and Figs. 4a - 4c show, for a 
given value of Am^ from the interval 3.0 x 10~^ eV^ ^ Am^ ^ 5.0 x 10~^ eV^ of interest the 
corresponding deformed pp z/g and recoil-electron spectra change very little with the change 
of sin^ 29 within the interval 10^'^ ^ sin^ 29 ^ 2x 10~^. Therefore Fig. 3a actually illustrates 
the pp spectrum deformations which are caused by the MSW transitions when 

fe has a value in the wider interval 0.4 ^ fe ^ 3.0. If Am^ 6 x 10"^ eV^ the pp neutrino 
flux suppression factor reads P^g^(E) = 1 — l/2sin^26', E < 0.42 MeV, and its deviation 
from 1 is unobservable. 

The Ue ^ i^s NA solution at 0.6 < f b ^ 2.0, for which 5 x lO^^ ^ sin2 2^ ^ 2 x 10"^ 
and 3.3 x 10~^ eV^ ^ Am^ ^ 8.0 x 10~^ eV^, implies (for given Am^ and sin^26') the same 
distinctive deformation of the *B z/g spectrum as the analogous z^g ^fj-ir) NA solution. 
Detailed studies of the *B z/g spectrum deformation due to the MSW transitions in the case 
of interest can be found in refs. [19,23,32] (see also Fig. 2a). However, the distortion of the 
spectrum of the recoil electron from the reaction u + e~ —>■ u + e' caused by neutrinos 
to be measured by the Super-Kamiokande collaboration is predicted to be larger than that 



suppression of the signal in BOREXINO due to the 0.862 MeV ^Be neutrinos by a factor 
not exceeding 0.15 concomitant with a strong distortion of the pp z^g spectrum is predicted in the 
case of the "low" $Bi 0.35 ^ f b ^ 0.43, vacuum z^g <^ Ug oscillation solution of the solar neutrino 
problem [22]. However, the z^g Vg oscillations would deform the lower energy part of the pp 
Ue spectrum, E ^ 0.30 MeV, and the distortions will be very different from those caused by the 
Ve MSW transitions (compare our Fig. 3a with Fig. 3b in ref. [22]). 
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implied by the similar u^. — > solution for the same Am^ and sin^ 29 (compare our Fig. 
2b with Figs. 6c, 6d in ref. [19] and with Fig. 3 in ref. [23]). By measuring the two 
spectra (in SNO and Super-Kamiokande experiments) and by testing the different types of 
correlations between them the and the Ve ^iJt{T) transitions lead to, it can be 

possible to distinguish between the MSW nonadiabatic Vg and v^, ^^(^t-) solutions 

for 0.6 < f B ^ 2.0. Other possible tests of the *B neutrino —>■ Vg conversion hypothesis, 
which could be performed using the SNO and Super-Kamiokande data provided the relevant 
conversion probability exhibits a nontrivial neutrino energy dependence (as in the case of 
the NA solution under discussion), were proposed in ref. [33]. Further, if 0.6 < f b ^ 2.0, 
the 0.862 MeV ^Be v^, flux is predicted to be suppressed due to the NA — > Vg transitions 
by a factor PmI^ < 0.15 for most of the values of Am^ and sin^ 1Q from the allowed region 
corresponding to given fB and fBe- For instance, for sin^ 2^ = 5 x 10^'^; 10~^; 2 x 10^^ 
we have P|^|^ < 0.15 if respectively Am^ ^ 7.3; 7.0; 6.8 x 10"^ eV^. For the solution 
under discussion P^gw takes values in the interval ~ (0.003 — 0.40), the maximal value being 
reached for Am^ = 8.0 x lO"*^ eV^ and sin^ 26* ^ 5 x lO'^. If Am^ > 5.0 x 10"^ eV^ the pp 
flux will be practically unaffected by the NA — > Vg transitions as Pmsw(E) — 1 — 1/2 sin^ 2^, 
while for smaller values of Am^ they will cause a rather large (and, perhaps, observable) 
distortion of the pp z/g and the corresponding recoil-electron spectra (Figs. 3a, 4a - 4d). 
For "large" values of $b, > 2.0, there exists a i/g — > NA solution, but not a similar 
— ^ix{T) one. The — > Vg solution holds for 10~^ ^ sin^ 2^ ^ 6 x 10~^ (Figs, la - 
Ic). Its distinctive features are: i) characteristic distortions of the ®B v^. (E > 6.44 MeV) 
and of the corresponding recoil-electron (Ee > 5.0 MeV) spectra (Figs. 2a and 2b), ii) 
Rcc/NC = R|c/NC> iii) Pmsw < 0.15 occurring for, e.g., sin2 2^ = 10"^; 2 x 10-^; 4 x 10"^ 
in the region Am^ ^ 7.0; 6.8; 6.7 x 10^^ cV^, and iv) substantial distortion of the pp Vf. 
(and the recoil-e~) spectrum for Am^ 4.5 x 10~^ eV^. Actually, with the change of fs, 
Am^ and sin^ 2Q within their respective intervals the 0.862 MeV ^Be neutrino suppression 
factor P|J^ is predicted to change from ~ 3 x lO"-'^ (at Am^ ^ 3.3 x 10"^ eV^) to ~ 0.34 
(for Am^ = 8.0 x 10~^ eV^). The minimal value of Pmsw depends strongly on sin^2^ 
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increasing from 3 x 10~^ to approximately 10~^ (6 x 10~^) when sin^ 26 changes from 10~^ 
to 4 X 10^^ (2 X 10^^). At the same time the maximal value varies insignificantly (decreasing 
by ~ 0.02) with the variation of sin^ 26. 

For the Ve ~^ AD solution, which we have found to exist at fe ^ 3.0 (which is strongly 
incompatible with the current solar model predictions [4-9] and which we shall consider 
briefly for completeness), the z/g flux will be suppressed by the energy independent factor 
P^g^ = sin^ 6 and there will be no distortions of the z/g and of the corresponding recoil- 
electron spectra. One will have also Rcc/nc — R-cc/nc- -'■^ contrast, in the case of the 
Ve v^(r) AD solution possible if 1.0 ^ fe ^ 3.4 [23], one has [19] Rcc/nc < 0.4 Rcc/nc B 
and the spectra of the ^B v^, having E > 6.44 MeV and of the recoil e~ with Ee > 5.0 MeV 
can be strongly deformed (see the first article quoted in ref. [18] and ref. [19]). The indicated 
strong spectra deformations are predicted to take place in the Ve and e~ energy intervals 
(6.44 - 11.0) MeV and (5.0 - 11.0) MeV, respectively, for Am^ ^ (4.5 - 5.0) x 10"^ eV^ 
and are clearly distinguishable [19] from, e.g., those associated with the ^ ^^t(T) NA 
solution at 0.6 < fe ^ 2.0 (see Figs. 4c, 5c and 6c in ref. [19]). For the z/g — > Vs AD 
solution the allowed values of Am^ and sin^ 2^ for, e.g., 3.0 ^ fe ^ 5.0 lie in the region 
5.0 X 10^^ eV2 ^ Am2 ^ 8.6 x 10"^ eV^, 0.26 ^ fiin^26 ^ 0.44. In this case we have 
0.15 ^ Pmsw ^ 0.37. The pp Ue flux will be suppressed by a factor (0.7 - 0.8), and the 
pp Ve and the associated recoil-electron spectra will be deformed for E > 0.22 MeV and 
To ^ 0.1 MeV (Figs. 3b and 5). Note that since the total pp neutrino flux $pp is tightly 
constrained by the data on the solar luminosity (see, e.g., [11]), an observation of a value of 
•I'pp < 0.85<l'p^, where <l>p^ is the pp neutrino flux predicted, e.g., in ref. [4], will be a very 
strong indication for an "unconventional" behavior of the solar electron neutrinos on their 
way to the Earth. 



^For values of Am^ ^ 2.0 x 10 ^ eV^ from the region of the Ve ^m(t) AD solution we have: 
Rcc/nc = sin^ 6 R^c/nc ^^^^ ^-^^ ^ 

sin 9 ^ 0.34. 
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There exist cosmological constraints on the MSW Ve — transitions [34] (see also [35]), 
which are relevant to the solution of the solar neutrino problem under discussion. They follow 
from the upper limit on the (effective) number of neutrino species, N^^, implied by big bang 
nucleosynthesis and hehum abundance considerations [36]. A limit of N^^ < 3.1; 3.4; 3.8 
is incompatible with the MSW adiabatic Vf. — > i/g transition solution for fe ^ 30; 9; 3. The 
cosmological upper limit on N^^ is presently a subject of debate [36] with values ranging 
from less than 3 to approximately 4. 

4. CONCLUSIONS 

The MSW i^s transition solution of the solar neutrino problem is stable with respect 
to changes in the predictions for the and ^Be neutrino fluxes: for 0.7 ^ fee ;^ 1.3 the data 
(1) - (4) admits nonadiabatic and adiabatic Vg — J^s transition solution for 0.4 ^ fe ^ 36 
and 3 ^ fe ^ 36 respectively. For 0.6 < fe ^ 36 the NA solution implies strong distortions of 
both the ^B i/g spectrum, measurable, e.g. in the CC disintegration of the deuteron in SNO, 
and of the spectrum of the recoil electron from the reaction i/ + e" — > u + e" induced by 
^B neutrinos, measurable both in SNO and Super-Kamiokandc experiments. This solution 
can be distinguished from the 1/^ — > NA solution at 0.6 < f b ;^ 2 (which implies weaker 
distortion of the spectrum of recoil electrons) by comparing the two spectra. Thus, if these 
spectra will be measured with the planned accuracy in the two experiments, they can be 
used not only to test the two solutions, but also to discriminate between them. 

Additional possible signatures of the — > Us NA solution under discussion are: i) a strong 
suppression of the signal due to the 0.862 MeV ^Be neutrinos in the BOREXINO detector by 
a factor PmIw ^ 0.15 (not possible in the case of, e.g., ^ti{T) MSW transitions), and/or 
ii) a specific and rather strong distortion (typically of the higher energy part, E ^ 0.30 MeV) 
of the pp Ve spectrum at E ^ 0.218 MeV and of the corresponding recoil-electron spectrum 
at 0.10 MeV ^ Te ;^ 0.26 MeV. If 0.4 ;^ fe ^ 0.6 these two signatures can be the only ones 
by which the — > Ug solution can be distinguished from the analogous Vf, — > MSW 
solution as in both cases the distortions of the ^B and the recoil-electron (Te ^ 5.0 MeV) 
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spectra can be unobservably small. Finally, the i/e — > t^s AD solution for 3.0 ^ fe ^ 5.0 
(which is a totally unrealistic possibility from the point of view of the solar models currently 
discussed in the literature) implies: i) an energy-independent reduction of the *B i/g flux by 
a factor 0.07 ^ sin^^ ^ 0.13, ii) a suppression of the BOREXINO signal due to the ^Be 
neutrinos by a factor (0.33 - 0.50), iii) a relatively strong suppression (by a factor (0.7 - 
0.8) of the total pp Ve flux, and iv) a characteristic but not strong distortion of the pp 
(recoil-electron) spectrum (Figs. 3b and 5). The last two signatures can be crucial for the 
identification of the — > Vg AD solution. 
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Figure Captions 

Figs, la — Ic. Regions of values of the parameters Aw? and sin^ 20 for which the solar 
neutrino data can be described at 95% C.L. in terms of Vg MSW transitions of the 

solar Vf, for values of fee equal to: a) 0.7; b) 1.0; c) 1.3, and for values of fe from the interval 
(0.40 - 35). Shown in the figures is also the region excluded at 90% C.L. by the absence of 
a difference between the day and the night event rates in the Kamiokande detector [12] (the 
region inside the dashed line) . The plot for this excluded region was taken from ref . [37] . 

Fig. 2a. Deformations of the spectrum of (i/g) neutrinos in the case of MSW transi- 
tions i^e — > fs of solar neutrinos for values of Am^ and sin^ 26 indicated in the figure. Each 
of the predicted z/g spectrum is divided by the standard one and the ratio so obtained 
is normalized to the value it has at the energy of i/g E = 10 MeV (see, e.g., ref. [19]). The 
error-bars shown illustrate the envisaged sensitivity of the data from SNO (CC reaction) 
after five years of operation; they do not include systematic errors. 

Fig. 2b. Deformations of the spectrum of e~ from the reaction v -\- v + e~ 

caused by ^B neutrinos, in the cases of MSW transitions Vg and for kinetic energy of 

the electron Te > 5 MeV. Each of the predicted recoil-electron spectrum is divided by the 
standard one and the ratio so obtained is normalized to the value it has at Tg = 10 MeV. The 
error-bars shown illustrate the sensitivity of the Kamiokande experiment [12] (large crosses) 
and the expected (ideal) sensitivity of the Super-Kamiokande experiment after five years of 
operation (small crosses); in the latter case no systematic errors have been included. 

Fig. 3a. Deformations of the spectrum of pp (t^e) neutrinos at E > 0.15 MeV in the 
case of the MSW Ve — ^ t^s "low" fe solution (see the text) for sin^ 2^ = 5 x 10~^ and several 
values of Am^. Each of the predicted pp v^. spectrum is divided by the standard one and 
the ratio so obtained is normalized to the value it has at the energy of the E = 0.3 MeV. 
The spectra shown in the figure change very little when sin^ 2^ is varied in the interval 
10-3 ^ sm^2e ^ 2 X 10-2 (see Figs. 4a - 4d). 

Fig. 3b. Deformations of the spectrum of pp {ve) neutrinos at E > 0.15 MeV in the 
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case of the MSW Ve — > fg adiabatic solution (see the text). The spectra are normahzed in 
the same way as those shown in Fig. 3a. 

Figs. 4a — 4d. Deformations of the spectrum of e~ from the reaction v -\- e~ ^ v -\- e~ 
caused by pp neutrinos, in the case of MSW Vf. Vs transitions (nonadiabatic solution) 
for sm^2e = 10'^ (a), 5 x lO'^ (6), lO'^ (c),2 x lO'^ (d), at Te > 0.10 MeV. Each of the 
predicted recoil-electron spectra is divided by the standard one and the ratio so obtained is 
normalized to the value it has at the kinetic energy of the electron Tg = 0.20 MeV. Note 
that the recoil-electron spectra change little with the change of sin^ 26. 

Fig. 5. Deformations of the recoil-electron spectrum due to the MSW Ug Vg transitions 
(adiabatic solution) at Te > 0.10 MeV and for sin^ 2d = 0.30; 0.40 and Am^ = 6.0; 7.0; 8.0 x 
10~^ eV^. The normahzation of the spectra is the same as in Figs. 4a - 4d. In each of the 
three pairs of sohd, dotted and dashed curves the upper (lower) curve at Te < 0.20 MeV 
(Te > 0.20 MeV) corresponds to sin^ 29 = 0.4. 



23 



